A glass anomaly in the shear modulus of solid "^He 
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The shear modulus of solid *He exhibits an anomalous change of order 10% [l|, at low tem- 
peratures that is qualitatively similar to the much smaller frequency change in torsional oscillator 
experiments. We propose that in sohd *He the stiffening of the shear modulus with decreasing 
temperature can be described with a glass susceptibility assuming a temperature dependent relax- 
ation time t(T). The glass susceptibility captures the freezing out of glassy degrees of freedom 
below a characteristic crossover temperature Tx ■ There the dynamic response of the solid satisfies 
uot{Tx) ~ 1, thus leading to an increase in the shear modulus. Within this model we predict that 
the maximum change of the amplitude of the shear modulus and the height of the dissipation peak 
are independent of the applied frequency uj. Our calculations also show a qualitative difference in 
behavior of the shear modulus depending on the temperature behavior of the glass relaxation time 
t(T). These predictions can be tested by comparing the complex shear modulus with experiments 
at different frequencies and for different levels of disorder. 
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The low-temperature anomaly of solid helium in tor- 
sional oscillators (TO) reported by Kim and Chan ^ has 
inspired an intense search for mechanical and thermody- 
namic anomalous properties. The observed increase of 
the resonance frequency below 200 mK was taken as 
evidence for the mass decoupling due to supersolidity - 
a quantum solid that can sustain mass superfiow with- 
out dissipation. It is now accepted that the presence of 
defects in solid ^He is required to produce supersolid sig- 
natures. Also it has been speculated that supersolidity 
may occur along dislocation lines [H| or grain bound- 
aries |7l49| in solid helium. Direct experimental evidence 
for a true phase transition into a supersolid state remains 
inconclusive. To date no definitive sign of Bose-Einstein 
condensation (BEC) has been seen in measurements of 
the mass flow H, H, [loll the melting curve [ll|, and the 
lattice structure [l^, ■ 

The basic issue with the unambiguous identification 
of the BEC features arises from the presence of defects 
that display their own dynamics and contribute to ob- 
servables in the same temperature and pressure range, 
where supersolid anomalies are expected. We therefore 
are left with the dilemma: on one hand defects are re- 
quired to produce supersolidity, on the other hand de- 
fects exhibit their own dynamics. Thus, any unambigu- 
ous identification of a possible supersolid state relies on 
a detailed understanding of the behavior of defects. For 
that reason we consider a theoretical framework that cap- 
tures the dynamics of defects in solid *He in the form 
of a glassy component, which makes up a small frac- 
tion of the crystal [iBl- This glassy component is 
suggested to cause the TO and thermodynamic anoma- 
lies. Further, it is consistent with reported signatures of 
long equilibration times, hysteresis, and a strong depen- 
dence on growth history. The glass may be created from 
distributions of crystal defects forming two-level-systems 
(TLS) [H, To determine the detailed nature of the 



TLS a detailed microscopic characterization of samples 
is needed. Possible candidates for the microscopic re- 
alization of the TLS are pinned segments of dislocation 
lines [H, [H, [ll] , which naturally occur in bulk and con- 
fined solid helium. These defects can be annealed away 
and hence drastically change the mechanical properties 
of the solid. In previous work [13, [IE fl8l - l20j we focused 
on the role of the glass on the TO and thermodynamic 
experiments. We demonstrated that the freezing out of 
defect dynamics can account for the observed anomalies. 
Further, we found that the predicted relaxation dynam- 
ics describes the observed low-temperature features when 
introducing a relaxation time t(T) that increases with 
decreasing temperature T. 



Here we consider the dynamic response of elastic prop- 




21|-|23j. Very recent shear 
reveal qualitative similar- 



erties in *He crystals 
modulus measurements ^ ^ 

ities with the TO experiments [1, [24- 30] . In the shear 
modulus experiment, the solid helium is grown in be- 
tween two closely spaced sound transducers. When one 
of the transducers applies an external strain, the other 
transducer measures the induced stress from which the 
shear modulus of the sample is deduced. The experi- 
ment thus provides a direct measurement of the elastic 
response to the applied force within a broad and tun- 
able frequency range. The frequency dependence of these 
elastic experiments can play an important role in deter- 
mining the origin of the observed anomalies. 

In this Letter we analyze the shear modulus within 
the same glass framework as was introduced for the TO 
and specific heat experiments ^13, [IB ll"8l - [20j . Within 
the glass model the amplitude of the shear modulus in- 
creases (stiffens) upon lowering T, because defects are 
freezing out. This is accompanied by a prominent dissi- 
pation peak, indicative of anelastic relaxation processes. 
These anomalies happen in the same temperature range 
where the TO anomalies were found. By studying the 
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glass model for a complex shear modulus ii{uj;T) we find: 
(1) The damping and amplitude of vibrations in "^He are 
controlled by the freeze-out dynamics of defects. They 
occur at temperatures where ujt{T) ^1. In our pic- 
ture the glassy contribution represents a small fraction 
of the total response because the glass occupies only a 
small fraction of the solid. Thus we propose a theoretical 
description of an elastic material with a small anelas- 
tic component that is modeled by a glassy susceptibil- 
ity. (2) We calculate the amplitude of the shear modu- 
lus, FIG. [TJ in agreement with shear mode experiments. 
Within this approach we find that the maximum of the 
shear modulus change and the height of the dissipation 
peak are independent of frequency. (3) We investigate 
both thermal (Vogel-Fulcher-Tamman) and nonthermal 
(power-law) activation processes. (4) We predict for both 
processes a relation for the inverse crossover temperature 
l/Tx vs. the applied frequency uj, i.e., ujt{Tx) = 1. 

Model - We investigate the dynamics of the solid jsil] in 
the presence of a back action following the logic laid out 
in Ref. 15|. Using standard response theory for elastic 



deformations we write the equation of motion 
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where p is the mass density and Ui is the displacement in 
the ith direction, ff'^'^ and ff^ are the external force 
density and the back action force density in the ith direc- 
tion, is the elastic stress tensor due to solid helium. 
In general. cr^° — Xijki Uki, with the elastic modulus ten- 
sor Xijki [3l|. The back action describes the delayed 
restoring force of a glass component that back-acts on 
the solid matrix and thus modifies the net force. For 
simplicity, we consider a homogeneous solid and set the 
shear wave propagation along the z axis and assume that 
the wave polarization lies in the x-y plane. For such a 
case the back action is 



BA 



dt'g{t,t'-T)dlu,{t') 



(2) 



where Q describes the strength of the back action on 
solid '^He and i = x,y. Although ff^ is typically 
much smaller than the elastic restoring force 9j(T^°, it 
is in fact this term that is responsible for the anomaly. 
The isotropic approximation is appropriate for measured 
polycrystalline and amorphous materials. 

Applying the same approximation to the elastic modu- 
lus tensor, Xijki = XoStjSki+poiStkSji+SuSjk), the elastic 
stress tensor in Eq. ([T]) is finite only for orientations j = z 
and either k or / equal to z. With k, I being interchange- 
able, the relevant element will be Xiziz, which gives fiQ. 
The fully dressed shear modulus relates the displacement 
to an external force, or pdfUi + pd^Ui,— ff''^^- Com- 
paring this expression with Eq. ([T]) we obtain 



p{uj;T) = po{T) - g{uj;T). 



(3) 



The back action may be described by a distribution of 
Debye relaxors G = dt P{t) [1 — 1/(1 — iujTt)] , with 



the dimensionless parameter t, the normalized distribu- 
tion of relaxation times P{t), and the relaxation time of 
the glass r. In general t{T) increases with decreasing T 
and approaches infinity at the ideal glass temperature Tg. 
The specific form of t{T) can change qualitatively the T- 
dependence of p and will be discussed. For simplicity, we 
choose for Q the Cole-Cole distribution. Integrating over 
Pit) yields 19]: 

g{LO-T) = go/(l-(*c.r)"), (4) 
M^;T) = p^[l~g/{l-{iLOTr)], (5) 

with the renormalized parameter g = go/poy which is 
sample dependent. The experimental measurables are 
the amplitude of the shear modulus, \p\, and the phase 
delay between the input and read-out signal, = arg (p); 
(j) measures the dissipation of the system, which is related 
to the inverse of the quality factor Q"^ = tan0. 

Several interesting results follow from Eq. (|4]). First, 
the change in shear modulus A/i between zero and infinite 
relaxation time is Ap/pQ = g, measuring the strength 
of the back action as well as the concentration of the 
TLS. Second, the peak height At/) of the phase angle is 
proportional to g. When ujt = 1 (location of dissipation 
peak) A0 becomes approximately 



g 



sin(Q!7r/2) 



(2 - .g)(l -cos(a7r/2)) 



(6) 



For 5^1 and 1 < a < 2, Eq. ^ simplifies even further: 
A(j) oc {2 — a)g(x (2 — a)Ap/po. The peak height A0 de- 
pends only on the phenomenological glass parameters a 
and g. So both Ap and A(j) are independent of frequency. 

Finally, we comment that a glass exhibits viscoelastic 
properties [32i | . So far the response formulation is equiv- 
alent to a viscoelastic system with a Cole- Cole distribu- 
tion of relaxation processes also known as a generalized 
Maxwell model ~ parallel connections of an infinite set of 
Debye relaxors each with a different r. 

Results - We now compare our glass model calculations 
with the experimental shear modulus measurements by 
Day and Beamish [l] for a transducer driven at 2000 Hz 
(SM2000), 200 Hz (SM200) and 20 Hz (SM20) 

To model specific examples of the T-dependent t{T) 
in Eq. ([4]), we consider: (1) the Vogel-Fulcher-Tammann 
(VFT) relaxation, and (2) the power-law (PL) relaxation, 
which represent thermal and nonthermal activation pro- 
cesses, respectively. For case (1) we assume the form: 



r(T) = 



To e 

CXD 



A/(T-T,) forr>r„ 



for T < Tg. 



(7) 



Here To is the attempt time and A is the activation en- 
ergy. The experimental data for all three frequencies are 
described by a single set of model parameters shown in 
FIG. [TJ Note that in our parameter search we did not 
bias Tg to be positive. It can be either positive (real 
transition at Tg, where r diverges) or negative Tg (no 
real transition). Good agreement between calculations 
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FIG. 1: (Color online) Experimental data and calculations 
(VFT) of the shear modulus vs. temperature. The red 
squares are the experimental data for the shear modulus am- 
plitude. The black-dash line shows the calculation of the 
shear modulus amplitude. The blue-solid line is the predic- 
tion for the phase angle (j). The calculation uses a — 1.35, 
g = 6.41 X 10"^^o = 16.2 fA Hz-\ tq = 39.3 ns, A = 1.09 
K, and Tg = -29.8 mK. The f=2000 Hz dataset is shifted 
by 0.13 fA Hz"^ to be comparable to the 20 Hz and 200 Hz 
datasets. 
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FIG. 2: (Color online) Experimental data and calculations 
(PL) of the shear modulus vs. temperature. The red squares 
are the experimental data of the shear modulus amplitude. 
The black-dash line shows the calculation of the shear mod- 
ulus amplitude. The blue-solid line shows the predicted 
phase delay. The model parameters are a — 1.45, g = 
6.32 X 10-^Ato = 16.2 fA Hz-\ to = 50.0 ns, p = 3.57, 
and To = 1.36 K. The f=2000 Hz dataset is shifted by 0.13 
fA Hz~^ to be comparable to the 20 Hz and 200 Hz datasets. 



and experiments are obtained for all three frequencies for 
Tg = -30 mK. We refer to this calculation as "VFT<". 

The calculated crossover temperature Tx is slightly 
higher than in experiment for the SM2000 and SM200 
datasets, while slightly below the SM20 dataset. In ad- 
dition, our model predicts the relative phase angle be- 
tween the input and output signal, which is related to 
the dissipation. The phase angle exhibits a peak at Tx, 
where the amplitude of the shear modulus has a turn- 
ing point. As expected Tx decreases with decreasing ui. 
By setting Tg = K the VFT expression reduces to an 
Arrhenius rule "VFTg". We find no distinct differences 
between the Arrhenius and non- Arrhenius VFT behavior 
(not shown). 

Next we examine case (2) of a power law for r: 

r(T) = TO {To/Tf, (8) 

which corresponds to a freeze-out temperature at T = 
K. To here is a parameter of temperature scale. The 
corresponding results are shown in FIG. [21 which we 
refer to as "PLq". . The calculation is almost equally 
good as the one labeled VFT< for datasets SM200 and 
SM20. However, for SM2000 the predicted crossover is 
broader and occurs at a temperature roughly 100 mK 
higher than in experiment. 

To search for a possible phase transition, we study 
the higher and lower frequency behavior for various re- 
laxation scenarios. For the Cole-Cole distribution the 
crossover happens when wt(T) = 1. From that we esti- 
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FIG. 3: (Color online). Prediction for the inverse crossover 
temperature vs. applied frequency. The black thick-solid 
(blue dashed) line is the VFT< (PLo) prediction, using the 
same parameters as in FIG. [T] (VFT) and FIG. [2] (PL). The 
black-thin line is the VFT prediction forcing Tg = K. The 
red dash-dotted line is the power-law prediction with a phase 
transition occurring at 40 mK (PL40), for which we used 
r = ro(|rg|/(r - Tg)Y for T > Tg and r = cx) for T < Tg. 



mate Tx as a function of the applied frequency / — uj/2n. 

Figure El shows 1/Tx vs. /. The VFT< calculation 
gives the best description to the experimental data. The 
VFT< line is convex, while the Arrhenius line VFTq is 
straight. The upward curvature is typical for a VFT 
relaxation time with Tg < 0, as opposed to Tg > 0. 
For comparison the power-law predictions are also shown 
for phase transitions occurring at K (PLq) and 40 
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FIG. 4: (Color online) The Cole-Cole plots for relaxation 
times VFT< and PLo (see text). For given form of r, all dif- 
ferent frequency curves collapse onto one single master curve 
reflecting that lot is the only scaling parameter. The VFT< 
curve is higher than the PLo curve because of its smaller expo- 
nent a. Both Cole-Cole plots show reflection symmetry about 
Re[|^ — /io|/A/i]=0.5, which is a consequence of the Cole-Cole 
distribution function. 

mK (PL40). For positive Tg (see PL40), we find a true 
freeze-out transition, indicating the arrested dynamics 
for / — > Hz. For both VFT and PL relaxation times 
our calculations demonstrate that in the low frequency 
limit the existence of a phase transition should show clear 
signatures of Tx converging toward the glass temperature 
Tg. This behavior can serve as experimental evidence for 
a possible phase transition. 

Finally, we show Cole-Cole plots for both relaxation 
times VFT< and PLq (see FIG. S]). In both cases, the 
calculated results for all three frequencies collapse onto 
a single master curve demonstrating that cut is the only 
scaling parameter in the glass model. Differences in mag- 
nitude between the VFT< and PLq calculations are due 
to different exponents a in the Cole-Cole distribution 
function of Debye relaxors. Of course, the use of other 
distribution functions may result in skewed Cole-Cole 



plots breaking the mirror symmetry. 

In summary, we have shown that the low-temperature 
shear modulus anomaly of solid "^He can be described 
using the theoretical framework of glasses. The elastic 
shear modulus is strongly affected by the dynamics of 
defects. The freezing out of defects leads to a stiffen- 
ing of the solid concomitant with a peak in dissipation. 
By studying the glass susceptibility due to the back ac- 
tion on solid helium, we find that both the amplitude 
change and T-dependence of the shear modulus are well 
captured by this model. An important consequence of 
the dynamic response analysis is the prediction of the 
dissipation or phase angle. In the proposed glass model, 
the peak height of the dissipation is independent of the 
applied frequency (in linear response) and linearly pro- 
portional to the Cole-Cole exponent a as well as the back 
action strength g. Since g depends on the concentration 
of the TLS, we predict that increasing disorder will result 
in larger amplitude changes of the shear modulus. Ad- 
ditionally, calculations of the inverse crossover temper- 
ature vs. applied frequency for different thermally acti- 
vated and power-law relaxation times show qualitatively 
different behavior. Finally, we find that for a positive 
glass temperature Tg the crossover temperature Tx con- 
verges toward Tg. This prediction can serve as a clear 
experimental demonstration for the existence of a true 
phase transition in solid *He at low temperatures. We 
hypothesize that the freezing out of fluctuating segments 
of dislocation lines are the relevant excitations contribut- 
ing to the reported anomalies in solid ^He. A detailed 
understanding of the microscopic nature of the glass and 
its excitations remains a pressing challenge. Future ex- 
perimental characterizations may elucidate the puzzles 
found in the dynamics and thermodynamics of solid ^He. 
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